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Abstract. We prove well-posedness of global solutions for a class of coagulation 
equations which exhibit the gelation phase transition. To this end, we solve an 
associated partial differential equation involving the generating functions before 
and after the phase transition. Applications include the classical Smoluchowski 
and Flory equations with multiplicative coagulation rate and the recently intro- 
duced symmetric model with limited aggregations. For the latter, we compute the 
limiting concentrations and we relate them to random graph models. 



1.1. Coagulation models. In this paper we deal with the problem of uniqueness 
of post-gelation solutions of several models of coagulation, namely Smoluchowski's 
and Flory's classical models, and the corresponding models with limited aggregations 
recently introduced by Bertoin [3|. 

Smoluchowski's coagulation equations describe the evolution of the concentrations 
of particles in a system where particles can perform pairwise coalescence, see e.g. 
[H [m [23]. In the original model of Smoluchowski [29], a pair of particles of mass, 
respectively, m and m', coalesce at rate K{m,m') and produce a particle of mass 
m + m'. In the discrete setting, the evolution of the concentration Ct{m) of particles 
of mass m G N* at time t > is given by the following system 

d 1 

—Ct{m) = -\ K(m, m') Ct{m') Ct{m — m') — \ K{m,m')ct{'m)ct{'m'). (1-1) 



Norris considered in [21] far more general models of cluster coagulation, where the 
rate of coalescence does not depend only on the mass of the particles but also on 
other parameters. In this general setting, most results on existence and uniqueness 
are obtained before a critical time, known as the gelation time, while the global 
behavior of the solutions after gelation, and in particular uniqueness, is not known. 

An example of a solvable cluster coagulation model is Bertoin's model with limited 
aggregations [3j, which we shall simply call the model with arms. In this case, 
particles have a mass but also carry a certain number of potential links, called 
arms. Two particles of mass m and m' may coagulate only if they have a positive 
number of arms, say a and a'. When they coagulate, an arm of each is used to 
create the bond and both arms are then deactivated, hence creating a particle with 
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a + a' — 2 arms and mass m + m'. The coagulation rate of these two particles is aa'. 
Therefore, if Ct{a,m) is the concentration of particles with a G N = {0, 1, . . . } arms 
and mass meN* = {l,2,...}, then the coagulation equation reads 



a+l m—1 



-Ct{a, m) = - a' [a + 2 — a')ct{a' , m')ct{a + 2 — a' ,m — m') 

a'=lm'=l 

— aa'ct{a,m)ct{a' ,m'). 



a'>l m'>l 



For monodisperse initial concentrations, i.e. co{a,m) = l{m=i}/x(a), with /i = 
(/^(ci))aeN a measure on N with unit mean, it is proved in [S] that this equation 
has a unique solution on some interval [0, T), where T = +oo if and only if K < 1, 
where 

K ■='^a{a-l)fi{a). (1.3) 

a>l 



In other words, if particles at time have, on average, few arms, equation (1.2) has 
a unique solution defined for all t > 0. When this is the time passes, all 

available arms are used to create bonds and only particles with no arms remain in 
the system. The limit concentrations Coo(0, m) as t — ?■ +oo of such particles turn out 
to be related to the distribution of the total population generated by a sub-critical 
Galton- Watson branching process (see e.g. 0) started from two ancestors: see [Hill] 
and section 11.41 below. 



1.2. The gelation phase transition. A formal computation shows that solutions 



of (1.1) with multiplicative kernel K{m,m') = mm' should have constant mass 

Mt ■.= ^mct{m), t>0, (1.4) 

m>l 

i.e. = 0. It is however well-known that if large particles can coagulate suf- 

ficiently fast, then one may observe in finite time a phenomenon called gelation, 
namely the formation of particles of infinite mass, the gel. These particles do not 
count in the computation of the mass so from the gelation time on, Mf starts to 
decrease. 



The reason why (1.2) can be solved, is that it can be transformed into a solvable 



PDE involving the generating function of {ct)t>o- In Equation (1.1), this transfor- 
mation is also possible for several particular choices of the kernel K{m,m'), namely 
when K is constant, additive or multiplicative: see e.g. 0. In the multiplicative case 
K{m, m') = mm', which is our main concern here, the total mass is a parameter of 
( 1.1[ ) and of the associated PDE, which is therefore easy to solve only when {Mt 



n)t>o 



is known. Existence and uniqueness of solutions of ( 1.1 ) are thus easy up to gelation, 
since in this regime, the total mass Mt is constant. 

After gelation, the gel may or may not interact with the other particles. If it 



does. Equation (1.1) has to be modified into Flory's equation (3.1). Else, the gel is 
inert, in which case Smoluchowski's equation continues to hold. Obviously, they are 
identical before gelation. 
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Occurrence of gelation depends heavily on the choice of the coagulation rate 
K{m,m'), and in the multiplicative case, gelation always occurs [TOl [121 HZ]- After 
gelation, the mass is not known, so Mt itself becomes an unknown of the equation, 
and well-posedness of the equation is then much less trivial. The multiplicative ker- 
nel is therefore particularly interesting, since it exhibits a non-trivial behavior but 
can still be studied in detail by means of explicit computations. 



The same phenomenon of gelation has been observed in [3] for (1.2) for monodis- 
perse initial concentrations cq. A formal computation shows that the the mean 
number of arms 

At '■= act{a,m), t>0, 

a,m>l 

satisfies the equation ^At = —A^ and should therefore be equal to for all t > 0. 
In fact, this explicit expression holds only until a critical time, which is shown to 



be equal to 1/{K — 1) ii K > 1 and to +oo if i^' < 1, where K is defined in (1.3). 
Again, the associated PDE is easy to solve before gelation since then. At is known, 
while afterwards, the PDE contains the unknown parameter At. 

1.3. Main result. In this paper we investigate the global behavior of Smoluchow- 



ski's equation with arms (1.2) before, at and after the gelation phase transition, 
proving existence and uniqueness of global solutions for a large class of initial con- 
ditions. The technique used, as in [3], is to transform the equation into a PDE. 
Since the total number of arms {At)t>{) is not a priori known, this PDE is non local, 
unlike the one obtained in the regime before gelation. This is the main difficulty we 
have to deal with. We use a modification of the classical method of characteristics 



to show uniqueness of solutions to this PDE, and hence to (1.2). We can consider 
initial conditions (co(a,m),a G N, m G N*) with an initial infinite number of arms, 
that is, such that 

Aq := aco{a, m) 

a,m>l 

is infinite, and show that there is a unique solution "coming down from infinity 
sufficiently fast", i.e. such that, for positive t, 

ft 

Al ds < +00. 



Note however that this is no technical condition, but a mere assumption to ensure 
that the equation is well-defined. 

We also consider a modification of this model which corresponds to Flory's equa- 
tion for the model with arms. In this setting, the infinite mass particles, that is, 
the gel, interact with the other particles. We also prove existence, uniqueness and 
study the behavior of the solutions for this model. 

In both cases, our technique provides a representation formula allowing to com- 
pute various quantities, as the mean number of arms in the system and the limiting 
concentrations. In Flory's case, we extend to all possible initial concentrations the 
computations done in [3] in absence of gelation. In the first model, a slight modifi- 



cation appears which calls for a probabilistic interpretation; see section 1.4 below 
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This seems to be the first case of a cluster coagulation model for which global 
well-posedness in presence of gelation can be proven. Another setting to which 
these techniques could be applied is the coagulation model with mating introduced 
in El. 



1.4. Limiting concentrations. In [3], explicit solutions to (1.2) are given for 
monodisperse initial conditions Co(a, m) = /i(a)l{m=i} for some measure /i on N 
with unit first moment. In particular, when there is no gelation, i.e. K < 1 where 



K is as in (1.3), and /i 7^ ^62, there are limiting concentrations 



m(m — 1) 



2)l|,=o}, m>2, 



where z/(m) = (m + l)/i(m + 1) is a probability measure on N different from 61. 
This formula clearly resembles the well-known formula of Dwass [7] , which provides 
the law of the total progeny T of a Galton- Watson process with reproduction law 
z/, started from two ancestors: 

P(r = m) = — z/*'"(m-2), m>2. 
m 

The similarity between the two formulas is no coincidence and is explained in [3] by 
means of the configuration model. For basics on Galton- Watson processes, see e.g. 

i- 

Let us briefly explain the result of [3], referring e.g. to [26J for more results 
on general random graphs. The configuration model aims at producing a random 
graph whose vertices have a prescribed degree. To this end, consider a number n 
of vertices, each being given independently a number of arms (that is, half-edges) 
distributed according to fi. Then, two arms in the system are chosen uniformly and 
independently, and form an edge between the corresponding vertices. This procedure 
is repeated until there are no more available arms. Hence, one arrives to a final state 
which can be described as a collection of random graphs. Then Corollary 2 in [4J and 
the discussion below show that, when there is no gelation, the proportion of trees of 
size m tends to Coo(0, m) when the number n of vertices tends to infinity. Hence, the 
final states in the configuration model and in Smoluchowski's equation with arms 
coincide. This shows that the former is a good discrete model for coagulation. 

Interestingly, the absence-of-gelation condition < 1 is equivalent to (sub)- 
criticality of the Galton- Watson branching process with reproduction law u, i.e. to 
almost sure extinction of the progeny, while K > 1 and gelation at finite time are 
equivalent to super-criticality of the GW process. 



In this paper, we obtain the limiting concentrations for (1.2) and its modified 
version when there is gelation. Let us start with the modified model, which is the 
counterpart of Flory's equation for the model with arms. In this case, and with the 
same notations as above, we obtain the limit concentrations 

Coo{a,m) = —-^ — — z/*™(m - 2) l{a=o}, m > 2, 
m{m — 1) 
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that is, the same exphcit form as the one obtained in absence of gelation. Again, this 
formula can be interpreted both in terms of a configuration model and of a super- 
critical Galton- Watson branching process. The relation between Flory's equation 
with arms and the configuration model is natural, since in both cases all particles 
interact with each other, no matter what their size is. It is worth noticing that, even 
though the limit concentrations have the same form with or without gelation, still 



some mass is eventually lost in presence of gelation, see (6.4) below. 

We also obtain the limiting concentrations for Smoluchowski's equation with arms, 
namely 

Coo(a, m) = ^ /3^-V*™(m - 2) l|,=o}, 
m[m — 1) 

where /3oo is some constant, which is 1 when there is no gelation, and is greater 



than 1 otherwise, see Section 6^ However, the probabilistic interpretation of /3oo is 
unclear. One can recover Smoluchowski's equation with arms from discrete models 
by preventing big particles from coagulating, as is done in [13] for the standard 
Smoluchowski equation, but the precise meaning of /3oo still seems to require some 
labor. 



1.5. Bibliographical comments. Smoluchowski's equation (1.1) has been exten- 
sively studied; we refer to the reviews [H [181 123]. Conditions on the kernel k are 
know for absence or presence of gelation, though this requires a precise definition of 
gelation, see e.g. [TT], or [13] in a probabilistic setting. For a general class of kernels 
Smoluchowski's solution has a unique solution before gelation [231 El [121 IB], and in 
the multiplicative case gelation always occurs [TUl Wl\ [T7]. 

For the monodisperse initial condition co(m) = l{m=i}, the first proof of existence 



and uniqueness to (1.1) before gelation is given in [20] , and a proof of global existence 
and uniqueness can be found in [15]. The case of general nonzero initial conditions 
has been considered by several papers in the Physics literature [8] [Q] [191 ISSl [31] , and 
by at least one mathematical paper |27J, which however treats in full details only 



the regime before gelation, see Remark 2.7 below. The same authors also provide 



in [2H] an exact formula for the post-gelation mass of (1.1), but with no rigorous 
proof. 



Thus, a clear statement about well-posedness of (1.1) for the most general initial 
conditions still seems to be missing, and our paper tries to fill this gap. We adapt 
the classical method of characteristics for generating functions, see [3 [3], which 
yields easily uniqueness before gelation for a multiplicative kernel [21]. We can in 
particular consider initial concentrations with infinite total mass, i.e. such that 

Mo := / mco(dm) = +oo, 

as long as /(o+oo)("^ ^ ^) co(dm) < +oo. This covers for instance initial conditions 
of the type Co(dm) = Cpm~'P dm with p G [1,2). 

Our main concern is uniqueness, since existence of solutions has been obtained in 
a much more general setting by analytic [Ml HTl 123] or probabilistic [131 E] means. 
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However, the case of an infinite initial mass seems to have been considered only in 



I6] in the discrete case, so we refer to Section 2.4 below for a proof. 



1.6. Plan of the article. We start off in Section 2 by considering existence, unique- 



ness and representation formulas for global solutions of (1.1), introducing and ex- 



ploiting all main techniques which are needed afterwards to tackle the same issues 



in the case of (1.2). We prove that for the most general initial conditions /io(dm), a 
positive measure on (0, +oo), Smoluchowski's equation with a multiplicative kernel 
has a unique solution before and after gelation. We also show existence and unique- 
ness for the modified version of Smoluchowski's model, namely Flory's equation, in 
Section 3. The techniques used are generalized in Section 4 and 5, where we prove 
analogous results for the models with arms. We compute the limiting concentrations 
in Section 6, which are not trivial, in comparison with the standard Smoluchowski 
and Flory cases, for which they are always zero. 

2. Smoluchowski's equation 



In this section we develop our method in the case of equation (1.1), proving 
existence, uniqueness and representation formulas for global solutions. Let us first 
fix some notations. 

• is the set of all non-negative finite measures on (0, +oo). 

• is the set of all non-negative Radon measures on (0, +oo). 

• For fi e and / G L^{^i) or / > 0, 



(/^,/) = / /("^) /i(dm). 

We will write m for the function m m, for m ^ m^, etc. 

For : (0, +oo) — )• M and m, m' > 0, A0(m, m') = (l){m+m') — (f){m) — (f){m'). 

Cc(0, +oo) is the space of continuous functions on (0, +oo) with compact 

support. 

For a function (t, x) h-)- (pt{x), (p[{x) is the partial derivative of with respect 
to X. 

d+ d+ 

— - or — — denotes the right partial derivative with respect to t. 
ot at 



We are interested in Smoluchowski's equation (1.1) with multiplicative coagulation 
kernel K{m, m!) = mm'. Note that the second requirement in the following definition 
is only present for the equation to make sense. 

Definition 2.1. Let fiQ G M.'^ ■ We say that a family {fit)t>o C A^^ solves Smolu- 
chowski's equation if 

• for every t > 0, (/is(dm), m)^ ds < +oo, 

• for all G Cc(0, +oo) and t > 

If' 

(/it, 0) = (/io, 0) + - / (/is(dm)/is(dm'),mm' A0(m,m')) ds, (2.1) 







2 

if (/io,wi^) < +00, then t h-)- (/xj,m^) is bounded in a right neighborhood of 
0. 
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The global behavior of this equation has been studied first for monodisperse initial 
conditions (i.e. /io = Si), in which case it can be proven that there is a unique 
solution (/it)t>o on IR+, which is also explicit, see [201 US]- This solution clearly 
exhibits the gelation phase transition. Up to the gelation time Tgei = 1, the total 
mass (/Ut,m) is constant and equal to 1, and then it decreases: {^t,fn) = lA 
t > 1. Moreover, the second moment (/it,m^) is finite before time 1, and then 
infinite on [l,-|-oo). It is also known in the literature that for any nonzero initial 
conditions, there is a gelation time < Tgei < +oo, such that there is a unique 
solution to (2.1) on [0,Tgei), and (/xt,m^) — +oo when t — )■ Tg"^: see e.g. 

Theorem 2.2. Let /io G Ai^ a non-null measure such that 



(/io,mAl)= / (m A 1) yUo(dm) < +00. (2.2) 

"'(0,+oo) 

We can then define 

Mo := (/io,m) G (0,+oo], K := {fio,m^) G (0,+oo], 
and the function 



go{x) := (/io,mx'^) = / mx™'/io(dm), x G [0,1] (2.3) 

J{0,+oo) 

with go{l) = Mo G (0,+oo]. Let 

Tgei:=l/iC G [0,+oo). (2.4) 



Then Smoluchowski's equation (2.1) has a unique solution on It has the follow- 
ing properties. 

(1) The total mass Mt = (/it,m) is continuous on [0, +oo). It is constant on 
[0,Tgei] and strictly decreasing on [Tgd, +oo). It is analytic on ]R+\{Tgei}. 

(2) If the following limit exists 

u:=- hm ^JfoM)!^e[-oo,0], 
x^i- gQ{x) +xgQ{x) 

then the right derivative M^^^^ of M at t = Tgei is equal to v. 

(3) Let rriQ = inf supp /io G [0, +oo). When t — )• +oo, 

1 



tM, 



(4) The second moment {^t,rn^) is finite for t G [0,Tgei) and infinite for t G 
[Tgci, +oo). 

Remark 2.3. • This result allows to recover the pre- and post-gelation for- 
mulas obtained with no rigorous proof in some earlier papers [21 El UHl [IHl 
[271 12H1 125] • The decrease of the mass in 1/t when mo > was also observed 
in these papers. Also, some upper bounds in 1/t for the mass were proven 

in [inin!. 
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• If m = 0, the mass tends to more slowly than small particles need to 
coagulate before any big particle can appear, and they coagulate really slowly. 
For instance, a straightforward computation shows that if /io(dm) = e~™dm, 
then Mt ~ More generally, the explicit formula in Proposition (2.6) 
allows to compute Mt for any initial conditions. 

• With this formula, it is easy to check that Mt^^^+ can be anything from —oo 
to 0. For instance, Mt^^i+ = for Qq^x) = {1 — x) log(l — x) + x, Mq = — oo 
for go{x) = a/1 — x log(l — x) + x, and for < a < +oo, Mq = —a for 
go{x) = 1 — ^/T^^x^. In particular, M need not be convex on [Tgei, +oo). 

2.1. Preliminaries. Let /iq be defined as in the previous statement. We shall 



prove that, starting from /iq, there is a unique solution to (2.1) on M"*", and give 



a representation formula for this solution. This allows to study the behavior of the 



moments. Let us start with some easy lemmas. So take a solution {fit)t>o to (2.1) 
and set 

Mt = {fit,m). (2.5) 

The two following lemmas are easy to prove, using monotone and dominated con- 
vergence. 

Lemma 2.4. {Mt)t>o is monotone non-increasing and right- continuous. Moreover, 
Mt < +00 for all t>0. 

Proof. Take (f)^{m) = m for m G [0,-?^], (f)^{m) = 2K — m for m G [K,2K], and 
(j)^{m) = for m > 2K, so that 0^ G Cc- Plugging (p^ in Smoluchowski's equation 



(2.1), letting K — +oo and using Fatou's lemma readily shows that {Mt)t>o is 
monotone non-increasing. Note also that t t— )■ Mt = sup {fit, (p^) is the supremum 
of a sequence of continuous functions and so is lower semi-continuous, which implies, 
for a monotone non-increasing function, right-continuity. Finiteness of Mt is now 



obvious since s i— )■ M^, and hence s h- )■ Mg, are integrable by Definition 2.1, □ 



Lemma 2.5. Assume that t i— )■ {fitjm'^) is bounded on some interval [0,Tq]. Then 
Mt = Mo forte [0,To]. 



By Lemma 2.4, {fit^m) < +oo for t > 0, so that we can define 



g^(x) = {fit,mx'^) = mx"" fitidm), xG[0,1], t>0, (2.6) 

"'(0,+oo) 

which is the generating function of mfit{dm). Then, using a standard approximation 
procedure, it is easy to see that g satisfies 

( gt{x) = go{x) + x{gs{x) - M,)^{x) ds, t > 0, x G (0, 1) 

I gtil) =Mt, t> 0. 

It is well-known, and will be proven again below, that Mt = Mq for all t < Tgei, since 
then, the PDF (2.7) can be solved by the method of characteristics: the function 



6i(x):[0,l]^[0,l] 



{X 



^^tiMo-goix)) ^ X e [0,1], t<T, 



gel 
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is one-to-one and onto, has an inverse ht : [0, 1] [0, 1] and we find 

gt{,x) = go{ht{x)), X e[0,l],t< Tgei. 

However Mt is not necessary constant for t > Tgei and tlie form of (pt has to be 
modified; we thus define 

= xate-*^°^''\ X G [0, 1], t > (2.8) 

where 

at := exp ( / M, ds ) , t > 0. (2.9) 



\"'0 / 

For t > Tgei, Mt is possibly less than Mq and 0j, which depends exphcitly on 
{Ms)s£[o,t], is possibly neither injective nor surjective. We shall prove that it is 
indeed possible to find it £ (0, 1) such that (f)t{x) : [0, 1] i— t- [0,it] is one-to-one and 
it is uniquely determined by Qq. 

2.2. Uniqueness of solutions. Using an adaptation of the method of characteris- 
tics, we are going to prove the following result. Note that in [27], this properties are 
claimed to be true but a proof seems to lack. We will use the same techniques in the 



proof of Theorem 4^ for the model with arms, but they are easier to understand in 
the present case. 



Proposition 2.6. Let {fit)t>o be a solution of Smoluchowski's equation (2.1). 

(1) For all t G [0,rgei], Mt = Mq = go{it), where it := 1. For all t > Tgei, 
Mt = go{it) where it G (0, 1) is uniquely defined by 

it9'o{it) = ]. (2.10) 

Moreover it and cj)t{-) satisfy 

<P't{it) = 0, <Ptm = 1 > U^), Vx G (0, 1). (2.11) 



(2) For all t > 0, the function <t>t{.') defined in (2.8) has a right inverse 

/ii:[0,l]^[0,£i], Uht{x))=x, XG[0,1], (2.12) 

and 

gt{x) = go{ht{x)), t>0, xg[0,1]. (2.13) 

(3) The functions {it)t>o md {Mt)t>o are continuous. 

(4) {fJ't)t>o is uniquely defined by fio- 

Remark 2.7. • For all t < Tgei, Mt = Mq, it = 1 and 0t : [0, 1] ^^ [0, 1] 
is one-to-one and onto. The first thing one needs to prove is that for all 
t > Tgei, it < 1, i-e. there is indeed x G [0, 1] such that 4>t{x) = 1, see Lemma 



2.9 the second one, is that it = mt, i.e. 0t(-) has an absolute maximum 



at it, see Lemma 2.10 In other words, one has to exclude the dotted lines 
as possible profiles of 0t(-) in Figure 1. These properties are not obvious, 
since (pt depends on {Ms)se[o,t] which is, at this point, unknown. All other 
properties are derived from these two. 
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1 t 



Figure 1. 0t before and after gelation. The dotted lines represent 
what (pt may look like. The solid one is the actual 0f. 

In [271 Section 6] one finds a discussion of post-gelation solutions, in par- 
ticular of the results of our Proposition 2.6 However this discussion falls 



short of a complete proof, since the two above-mentioned properties are not 
proven. In particular, no precise statement about what initial conditions can 
be considered is given. 

The following lemma is a list of obvious but useful properties satisfied by g and 



Lemma 2.8. The function g defined in (2.6) satisfies the following properties. 

(al) (t, x) ^ gt{x) is finite and continuous on [0, +oo) x [0,1); 

(a2) For all x E [0, 1), t gt{x) is right differentiable on (0, -|-oo); 

(a3) For all t > 0, x 9t{.x) is analytic on (0, 1) and monotone non- decreasing; 

(a4) For all t > 0, x gt{x) G [0, -|-cxo] is continuous on [0, 1]. 

The function defined in ( |2.8[ ) satisfies the following properties. 
(bl) (pt is continuous on [0,1] and analytic on (0,1); 
(b2) 0^(0) = 0, 0t(l) = e-/o(A^o-M.)ds ^ [0, 1]; 
(b3) (j){{x) = ate-*5«(^)(l - txg'o{x)) for x E (0, 1); 

(b4) For t < Tgei, (pt is increasing. For t > Tg^i, x xg'^lx) is increasing, 
(pt(0) > and 0^(1) < 0. In particular, for t > Tgei, there is precisely one 
point rrit G (0, 1) such that 

ci>'t{mt) = 0. (2.14) 

(b5) For t > Tgei, (pt is increasing on [0,mi] and decreasing on [rrit, 1]. 
Moreover, 

(cl) The map {t,x) i— )■ 0t(x) is continuous on IR+ x [0, 1); 
(c2) The map {t,x) h- i- 4>i{x) is continuous on IR+ x (0, 1); 
(c3) For every x G [0, 1), t 0t(x) is right differentiable and 

^ = 0t(x)(Mi - go{x)) X G [0, 1), t > 0. (2.15) 



POST-GELATION UNIQUENESS FOR COAGULATION EQUATIONS 



11 



Property (b5) implies that there are at most two points in (0, 1) where ipt equals 
1. Take it to be the smallest, if any, i.e. 

it = inf{x > : = 1} (inf := 1). (2.16) 

Lemma 2.9. (1) For every t >0 and every x G [0,it] 

gtiU^)) = go{x). (2.17) 

(2) For all t G [0, Tgd], It = 1? and for t > Tgd, < < 1. In particular, for all 
t>0, (f)t{it) = 1 and 

g,{it) = gt{l) = Mt. (2.18) 

(3) Finally, t it is monotone non-increasing and continuous on M+. 



Proof. (1) Let us first prove that there exists r > such that (2.17) holds for 
t E [0, r[. Fix < a < 6 < 1. Since < min[a_fe] (po < max[a,fe] (po < I, then by 
property (cl) there is r > such that 

< min (j)t < max0t < 1, Vt G [0,r). 

[a,b] [a,b] 

So, for a fixed x G [a, b], the function 



is well-defined and using (2.7) and (2.15), we see that 

«t=^ (^^|^(0s(a;)) + ^(0s(a;))^^(x)^ ds = j.Usds 
where 



Since x G [0, 1), supigjo^) |7t| < +oo and therefore Ut = 0. Hence (2.17) 
holds for X G [a,b] and t G [0,r[. Since both terms of (2.17) are analytic 
functions of x on {0,it), by analytic continuation, (2.17) actually holds on 
(0,^i), and hence on [0,it] by continuity. 
(2) Let us now extend this formula to t G M^. Let 

T = sup{t>0: VsG [0,t], Va; G [0,4], gsiM^)) = 9o{x)} > r > 0, 

assume T < +oo, and denote by £ the left limit of {it)t>o at T. First, i 
cannot be 0, since otherwise we would get when s — )■ T~ 

For every t < T~ , < £ < £t, so for every x G (0,£), gt{4>t{x)) = go{x) and 
4>t{x) < 1. Using the continuity property (cl) and passing to the limit when 
t — )• T~ in this equality, we get 

griM^)) = 9oix), Vxg(0,£). 

By the same reasoning as in point (i), we obtain a. T' > T such that 
gt{4>t{x)) = go{x) for all t G [T, T') and x in a non-empty open subset of 
(0,£). By analyticity and continuity, the formula gt{(f)t{x)) = go{x) holds for 
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every t G [T, T') and x G [0,£(]. This contradicts the definition of T, and so 
T = +00. This concludes the proof of point (1) of the Lemma. 
(3) For the statement (2) of the Lemma, let us show first that (/Xt, m?) is bounded 
on [0,To), for every Tq G [0,Tgei). Let T', the smallest time when this fails 
(provided of course th at Tgd > 0). By assumption (see Definition 2.1), 
T' > 0. Differentiating ( |2.17[ ) with respect to x and having x tend to £t = 1 
gives, for t < T', 



(5) 



1 



tK 



This quantity explodes only when t = Tgei = 1/K, so T' = Tgei. 
(4) The boundedness of {{fit,m'^))ti^io,To) just proven for all Tq G [0,Tgei) and 
2^ imply that for t G [0,Tgei), Mf = Mq. By the definition of 
lows that 0t(l) = 1 for t G [0,Tgei). But (pt is increasing, so £t = 1 
for t G [0,Tgei). Assume now that for some t > Tgei, = 1- Then (2.17) 



Lemma 

(pt, it fo: 



holds on [0,1], and this is impossible since the right term is an increasing 
function of x, whereas the left one decreases in a left neighborhood of 1 since 
(p'til) < 0. The fact that (pti^t) = 1 follows then directly from the definition 

Finally, the inequality it > is obvious 



of it and the continuity of 
since 0t(O) = 0, and computing (2.17) at x 



gives (2.18). This concludes 



1 and Mt = Mq for all t < Tgei. Now, go is strictly 



the proof of (2). 
We know that it 

increasing and continuous. Since {Mt)t>o is monotone non-increasing and 
right-continuous by Lemma 2.4, so is {it)t>o by (2.18). To get left-continuity 
of (£t)4>Tgcp consider t > Tgei, and let i be the left limit of is at t. We have 



< 



t+(t-Tg^i)/2 < 1, SO by the continuity property (cl) above. 



Hence 4>t{i) = 1- Assume i > it (that is, i is the second point where 4>t 
reaches 1). Take x G {it,i)- By property (b5), 0t(x) > 1. But on the other 
hand, x < i < is for s < t, so 0s(x) < 1, and so (pt{x) < 1, and this is a 
contradiction. So i = it and {it)t>o is indeed continuous. This concludes the 
proof of (3) and of the Lemma. 

□ 



Finally, we will see that for t > Tgd, it = mt, so that (pt increases from to 
1, which is its maximum, and then decreases. To this end, recall that {it)t>o is 
monotone non-increasing and that {it)t>o and {(pt)t>o are continuous, so the chain 



rule for Stieltjes integrals and (2.15) give 



1 = <Pm = M^o) + f<P's{Q dis + J' 

1+ / 01(4) d4+ / (ps{is){Ms-go{is))ds 











POST-GELATION UNIQUENESS FOR COAGULATION EQUATIONS 



13 



that is, with p.lSh 



Hence, d^ra.e (p't{it) = 0, i.e. 



(j)[{et) d£t = 0. (2.19) 
= rrit. This is actually true for all t > Tgei, as 



namely part 3 of the proof of Theorem 4.2 



we shall now prove. This result also has its counterpart in the model with arms, 

0, i.e. it = rrit, the point where (j)t 

Vt > Tgei. (2.20) 



Lemma 2.10. For every t > Tgei, 0t(^t) 
attains its maximum. In particular, 

1 



t 



Proof. First, recall that 



is increasing on [0,it], so that (l)t{it) > 0, that is 
1 



t 



(2.21) 



Assume now that there is a t > Tgei such that (f)'t{it) > 0, and consider 

s = sup{r G (Tgei,t) : 0'^(£,) = 0}. 

As noted before, t H- is strictly decreasing for t > Tgei for any t > Tgei, so 
d£t([Tgei, Tgei + e[) > for all e > 0. Hence there are points r < t where 4>'ri^r) = 0, 
and thus the definition of s does make sense. 

Take now (r„) a sequence of points such that T < r„ < t, 4>r{ir) = and (r„) 
converges to s. Since < < 1, by property (c2) above, we get 

o = C(4j^0;(4) 

so that 4>'s{is) = 0. This shows that s <t, and that for r G (s,t), (p'ri^r) > 0. Hence, 
by continuity of {ir)r>o and by ( |2.19[ ), {£r)re[s,t] is constant. This gives 



< 



which is a contradiction since s < t. In particular, (p'tiit) = implies (2.20). □ 
Proof of Proposition 2.6. By Lemma 2.10, necessarily Mj = Mq on [0, Tgei] and for 



t > Tgei, Mt := gt{l) = goiit), where 



(2.22) 



Since x i— )■ xg^i^x) is strictly increasing from [0, 1] to [0,/^], where K = {fio,m'^) = 
1/Tgei, this equation has a unique solution for t > Tgei. Hence Mt is uniquely defined. 
Therefore at and (pt are uniquely determined by qq, so we can define (pt as in (2.8), 
and Lemma [2^ shows that gt{(t)t{x)) = go{x) for x G [0, it], and that cpt is a bijection 



from [0,it] to [0,1]. So it has a right inverse ht, and compounding by ht in the 
previous formula gives 

gt{x) = go{ht{x)) (2.23) 
for all X G [0, 1], t > 0. Thus gt can be expressed by a formula involving only go and 
in particular, {fit)t>o depends only on fiQ. This shows the uniqueness of a solution 
to Smoluchowski's equation (2.1). □ 
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2.3. Behavior of the moments. In this paragraph, we will study the behavior 
of the first and second moment of {^t)t>o as time passes, showing how to prove 
rigorously and recover the results of [9]. For more general coagulation rates, one can 
obtain upper bounds of the same nature, see [17]. 

First consider the mass Mt = {nt^rn). We will always assume that Tgd < +oo. 
Let us start with the following lemma. 

Lemma 2.11. Let v E Ai^ he a measure which integrates x ^ for small enough 
y > 0. Let mo be the infimum of its support. Then 

{iy,xy'') 
lim — — = mo. 

Proof. First, note that xy^ > my^ u-a.e. so 

{u,xy^) 
limmi — — > mo. 

y~^o {u, y^) 

Let us prove now that 

hmsup — < mo. 

y^o (z^, yn 

Assume this is not true. Then, up to extraction of a subsequence, we may assume 
that there exists a > such that for arbitrary small y G (0, 1), {y,xy^) > (mo + 
a){v,y^). Hence (z/, (x — mo — Ci)y^) > 0, so 

(i/, (x - mo - a)?/''l{^>„o+„}) > (i/, (mo + a - x)y''l{jno<x<mo+a}) ■ (2.24) 

But 

(Z/, (mo + a - x)y'^l{mo<x<mo+a}) > [mo + a- x)y''l{mo<x<mo+a/2}) 

> (Z/, (mo + a - x)l{mo<x<mo+a/2})y""'^°'^'^ 

and 

(Z/, (X - mo - a)y''l{^yrno+a}) < {v, {x-uiQ- a)l{x>mo+Q})y'"°"^"- 



With ( |2.24D , this shows that 

(Z/, (X - mo - a)l{^>rno+a})y''^'^ > (mo + « - x)l{mo<x<mo+a/2}) 

and having y tend to zero gives 

> (Z/, (mo + a - x)l{mo<x<mo+a/2}) 

which is a contradiction since z/([mo,mo + tt/2]) > 0. □ 

Corollary 2.12. The mass of the system is continuous and positive. It is strictly 
decreasing on [Tgci,+oo). Moreover, denote mo = inf supp /io- Then 

lim = mo. 

t^ + OO tMf 
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Proof. Recall that Mt = Qoi^t) so the first properties follow from Lemma 2.9 Denote 
now i/(dm) = m/io(dm). For t > Tgd, itg^^it) = 1/^, so 

1 {u,xi^) 



and since £j — )■ when t — )■ +oo, this tends to mo by Lemma 2.11 



□ 



We can also study the behavior of the mass for small times. Recall that before 
gelation, the mass is constant at 1. We have seen that it is continuous at the gelation 
time. We may then wonder if its derivative is continuous, that is if Mt^^^^ is zero or 
not. 

Lemma 2.13. The right derivative of M at Tgei is given by 



Mt 



provided the limit exists. 



G [-00,0] 



Proof of Lemma\2J^ For t > Tgei, f{it) = l/t with f{x) = xgQ{x), and < it <1. 
But f'{it) 7^ 0, so by the inverse mapping theorem, {£t)t>o is different iable and 



t'f'iitY 

Using the fact that Mf = goi^t), it is then easy to see that 



M, 



g',{lt)+Wt 



Since {(!-t)t>o is continuous at Tgei and £^^1 = 1; the result follows. 



□ 



Recall that the gelation time is precisely the first time when the second moment 
(/i4,m^) of {fit)t>o becomes infinite. It actually remains infinite afterwards. 

Corollary 2.14. For all t > Tgei, (fJ^tyrri^) = +oo. 

Proof. Note that 

(/if,m^) = g[{l), 



this formula being understood as a monotone limit. By (2.17), for x < it 

<P't{x)g't{(l)t{x)) = g'oi^)- 

When x^£t, by Lemma 2.10 and g'^^x) g'o{it) 7^ since it > 0. So 

g'ti<Ptm = g'tii) = +00. 



□ 
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2.4. Existence of solutions. Existence of solutions of (2.1 ) is a well-known topic, 
flBI . However, the case Mq 



see e.g. [idj. However, tne case ivIq = +00 is apparently new, so that we give a 
short proof for the general case based on previous paper s, m ainly p7]. 

Let now /iq G A^j be as in the statement of Theorem 2.2 and let us set Qq as in 
(2.3), it and Mt as in point (1) of Proposition 2.6, at and (ft as in (2.9) and (2.8). 
Then it is easy to see that (pt admits a right inverse ht satisfying (2.12), and we can 
thus define 

gt{x) := go{ht{x)), t>0, xE [0,1]. 

It is an easy but tedious task to check that gt satisfies (2.7) and all properties (al)- 
(a4) above. In particular, if go{l) = +00 then ht{l) < 1 and therefore gt{i) < +00 
for all t > 0. Following [27], we can now prove the following. 



Proposition 2.15. For all t > there exists fit E Ai^^ such that 



gt{x) = {fit,mx'^) 



mx"^ fit{dm), X e [0, 1]. 



'(0,+oo) 

Proof. Let t > be fixed. We set for all ?/ > 



^{y)-=tgo{e ^), G{y) ■.= T{y) + y -log at = -log (f)t{e~ 



We recall that / : [0, +00) t-)- [0, +00) is completely monotone if / is continuous on 
[0, +00), infinitely many times differentiable on (0, +00) and 

(-l)'0(y)>O, VA;>0,yG(0,+cx)). 

It is easy to see that $ and F are completely monotone. Moreover, G has a right 
inverse 

: [0,+oo) ^ [log(l/£i),+oo), G-\y) = - log /it(e-^), y>0, 
and therefore by the definitions 



goihtie-y)) = HG-\y)) 



y>o. 



By [271 Thm. 3.2], $ o G ^ is completely monotone and therefore, by Bernstein's 
Theorem, there exists a unique z/^ G A^j such that 

9t{e-y) = g,{ht{e-y)) = m-\y)) = [ e-^™z/,(dm), y>0. 

Since gti^) < +00 for all t > 0, we obtain that {ut^m) < +00, so that we can set 
Utidm) := mz/i(dm), and we have found that there is a unique Ht E M.^ such that 



9t{x) = goihtix)) 



x'^ m fit{dm) 



X G (0,1]. 



(0,+oo) 



□ 



In order to show that (/ii)t>o is a solution of Smoluchowski's equation in the sense 
we have to check that dt < +00 for all e > 0. This is the 



2.1 



of Definition 
content of the next result. 
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Lemma 2.16. // {iJ,t)t>o is the family constructed in Proposition 2.15, then for all 
e > 0, f^i^J's, ^)'^ ds < +00. 

Proof. If Mq < +00 then there is nothing to prove, since {Mt)t>Q is monotone non- 
increasing, so let us consider the case Mq = +00 and thus Tgei = 0. Since = go{it) 
is bounded and continuous for t G [6, e] for all 6 G (0, e), we have by (2.10) and (2.3) 



M^dt 



gl{£t)dt = egl{e,)-6gl{es 



< eg^oiQ - I 2 



egoi^e 



2tgo{it)g'oiit)dit 

9o{y) 



< egoiie) + -r (/io, 



y 

Tfi 2 

) < e^o(4) + Y (/io, m A 1). 



Letting 5 | 0, by (2.2) we obtain the desired result. 



□ 



We now finish the proof of existence of a solution by showing that {^t)t>o indeed 
solves (2.1). By choosing x = e~^, ?/ > 0, in (2.7), we find an equality between 
Laplace transforms. Since the Laplace transform is one-to-one, then we obtain 
(0. 

Remark 2.17. In the proof of uniqueness, we may only require that {fiQ,my"^) < 



+00 for every y G [0, 1). However, the same kind of computation as in Lemma 2.16 
shows that if this the case, but {fio,m A 1) = +00, then Jq M"^ ds = +00 for all 
t > 0, in contradiction with Definition 12. 11 of a solution. 



3. FLORY'S EQUATION 

We will now consider the modified version of Smoluchowski's equation, also known 
as Flory's equation, with a multiplicative kernel. 

Definition 3.1. Let /io G M.^- ^^V ^^^^ ^ family {fit)t>o C Ai^ solves Flory's 
equation (2.1) if 

• for every t > 0, (/^^(dm), m)^ ds < +00, 

• for all G Cc(0, +00) and t > 



+ 



{fis{dm)fis{dm'),mm' A(j){m,m')) ds 



(3.1) 



{fis,(p){fio{dm) - /i^(dm),m) ds. 



if {jjiQ^m?) < +00, then t 1— )■ {fj,t,m'^) is bounded in a right neighborhood of 
0. 



In equation (3.1), the mass that vanishes in the gel interacts with the other par- 
ticles. It is a modified Smoluchowski's equation, where a term has been added, 
representing the interaction of the particles of mass m with the gel, whose mass is 
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i.e. precisely the missing mass of the system. Notice that in this case the equation 
makes sense only if (/io,m) < +00. 

The mass is expected to decrease faster in this case than for (2.1 ). This is actually 
true, as we can see in the following result. 

Theorem 3.2. Let fiQ G Ai'^ a non-null measure such that {^q^itl) < +00, and set 
Mo := (/io,m) G (0,+oo), K := {fio,m^) e (0,+oo]. 



Let Tgei := 1/K G [0, +00). Then Flory's equation (3.1) has a unique solution 
(/it)t>o on ]R+. It has the following properties. 

(1) We have Mt = go{k), where k = 1 fort < Tgei and, fort > Tgei, k is uniquely 
defined by 

l^ = e-'{M,-9,(U))^ G[0,1). 

Therefore t ^ Mt is continuous on [0,+oo), constant on [0,Tgei], strictly 
decreasing on [Tgei, +00) and analytic on ]R+\{Tgoi}. 

(2) The function (pt{x) = xe*^*'^''"^"'-^-''' has a right inverse ht : [0,1] [0,/t]. 
The generating function gt of {^t)t>o is given for t > by 

9tix) = go{ht{x)). 

(3) Let mo = inf supp fiQ > 0. Then, when t — )■ +00, 

M,e"^°* ^ mo^o({mo}) 

and for every e > 

(4) The second moment {m'^,Ct) is finite on ]R+\{Tgci} and infinite at Tgei. 

Remark 3.3. • Norris Thm 2.8] has a proof of global uniqueness of 
Flory's equation (3.1) for slightly less general initial conditions (/io such that 
{fiQ, 1 + m) < +00), but for a much more general model. 
• When mo > 0, it was already observed (Proposition 5.3 in [^) that the 
mass decays (at least) exponentially fast (see also [HI ESI EI])- 



Proof of Theorem \3.^ The proof is very similar to (and actually easier than) that 
of Theorem 12.21 



(1) Arguing as in the proof of Lemma 2.4, we obtain easily that {Mt)t>o is 
monotone non-increasing and right-continuous. As in Lemma 2.5, if t i-> 
{fit, m?) is bounded on some interval [0, To], then Mt = Mq for t G [0, To] and 
therefore {fit)t>o is a solution of Smoluchowski's equation (2.1) on [0,To]. 

(2) Consider initial concentrations /xq as in the statement, a solution {fit)t>o to 
Flory's equation and gt{x), x G [0, 1], generating function of m jjit{dm). Then 
gt solves the PDE 

^ = x{gt-Mo)^, Vt>0,xG[0,l], (3.2) 
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the same as the one obtained for Smoluchowski's equation before gelation. 
It may be solved using the method of characteristics. Indeed, the mapping 

(j)t{x) = a;e*(^'°-^o(")) =x+ [ (Mq - go{x)) ds, (3.3) 

Jo 

has the following properties 
(dl) 0i(O) = 0, 0t(l) = 1. 

(d2) For all t > 0, (p[{x) = e*(*'^o-5oW)(l - txg'^ix)). 

(d3) For if: < Tgei, 0t(-) is increasing; therefore, 0t(x) G [0, 1] for all x G [0, 1] 

and 0t(x) = 1 if and only if x = 1 
(d4) For t > Tgei, 0t(-) is increasing on [0, mt] and decreasing on [m^,!], 

where nit is the unique x G (0, 1) such that 4>t{x) = 0, i.e. such that 

ixg^ix) = 1. 

(d5) For t > Tgei, (pt{^t) > 1, since 0t(l) = 1 and 4>t{l) < 0. Therefore there 

is a unique k G {0,mt) such that 4>t{k) = 1- 
(d6) For t > Tgci, (t>'t{lt) 7^ 0, since It < rrit. 















1 


t 


1 



Figure 2. (j)t before and after gelation. 
Setting It := 1 for t < Tgei, 0t is thus a continuous bijection from [0,/^] to 



(3.3) and arguing as in part (i) and (ii) of the proof of Lemma 2.9 

= 9Mt{x)) 



[0, 1], with continuous inverse function ht : [0, 1] (-)■ [0, lt\. By using (3.2) and 

roof of Lemn 

goi^x) satisfies Ut{x 



we can 
Mo(x) = 



see that the function ut{x) 
for all t > and x G [0,/t]. Therefore the only solution of the PDF (3.2) is 
given by 

gt{x) = g^{ht{x)), t>0, xG[0,l]. (3.4) 

Flory's equation has thus a unique solution on M"*", and its generating func- 
tion is gt- 

(3) We have seen in (d5) above that, for t > Tgei, there is a unique It G [0, 1) 
such that (ptih) = 1- The relation (ptilt) = 1 with It G [0, 1) is equivalent to 
It = e"*(^^o~^o('*)) with It G [0, 1). This relation implies that t k is analytic 
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for t > Tgci- A differentiation shows that 
dk (Mo - goilt))k 



< 0, t>T, 



geb 



dt 1 - tg'^{lt)lt 

since g'Q{k)h < g'o{^t)'m't = ^/t and go{k) < fi'o(l) = Mq. Let i be the hmit 

r'gclf 



of It as 1 1 Tgei: then we obtain i = e-^g-i(^^o-foW), i.e. (pr = 1- By (d3) 



above, this is equivalent to £ = 1. 

(4) Since Mt = gtil) = go{ht{l)) = go{lt), the properties of t Mt = go{lt) 
follow from those of t It- Recall now that (ptih) = 1, that is 

\og{k) = t{go{lt)-l). (3.5) 

If the limit Z of as t — > +oo were nonzero, then passing to the limit in this 
equality would give log(/) = — oo. So / = and 

log/t~-l (3.6) 

• Assume m > 0. Now, obviously go{x) < x™, so 

\og{tgo{lt)) = log/t + log5(o(/f) < \ogt + mlogk -oo. 

Hence tgoih) ^ and yields logk + t ^ 0. Hence /J" ~ e"™*. 
Finally 

lim Mte^' = lim ^ = m^o({m}) 

since by dominated convergence, go{x)x^™' — )■ m/io({'^}) when x — 0. 
Now, by monotone convergence, if m' > m, then gQ{x)x^"^ — )■ +oo 
when X tends to 0, whence 

lim M^e™'* = lim ^ = +oo. 

t-s>+oo t-s>+oo /™ 

• Assume now m = and let e > 0. By monotone convergence gQ{x)x''^ — )■ 
+00 as X 4 0, so using (3.6) we see that g{lt)e~^^ — )■ +oo as t f +oo, 
which is the desired result. 

(5) Finally, (3.4) gives for x < 1 and t > Tgei 

When X t 1, ht{x) ^ k < I, and (pt{ht{x)) (t)'St) ^ by (d6) above. So 
{Ht,nn?) = g'til) < +oo. 



(6) Existence of a solution of (3.1) follows arguing as in section \2A 



□ 

Corollary 3.4. Let fiQ G A^^ such that {jj,Q,m) < +oo and let ifJ'f)t>o and (/xf)t>o 
the solutions of (2.1), respectively, (3.1). Then 

• fif = for all t < Tgei := l/(/^o, "n^^); 

• (yuf , m) < (/if, m) /or a// 1 > Tgei. 
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Proof. For all t < Tgei, {fi^,m) = ,m) and therefore /if solves (2.1), so that by 
uniqueness of Smoluchowski's equation we have that /if = /if. For t > Tgei we have 
Qoilt) while (/if,m) = go{it), where It and it are defined respectively 

by 



that (/if, m) 



and 



-t(Mo-S()(/t)) 



G [0, 1) 



In points (d4) and (d5) of the proof of Theorem 3.2, we have shown that k < rrit 
where tmtgQ{mt) = 1, so that nit = it < k- Hence {fif,m) = go{it) < go{k) = 
{l^^,m). □ 

As anticipated, the mass decreases faster in Flory's case than for Smoluchowski's 
equation. In particular, in Flory's case {fit, rn"^) becomes finite immediately after 
gelation, the mass remaining however continuous (we can think that the big particles, 
which have the biggest influence on this second moment, disappear into the gel). 
Moreover, if inf supp /Iq > then the mass decays exponentially fast, which is to be 
compared with the slow decrease in 1/t in Smoluchowski's equation. 

Remark 3.5. The mass in Flory's equation may decrease slower if inf supp /io = 0. 
For instance, if /io(dm) = e~''^dm, then Mt ~ t~^. 

4. The model with limited aggregation 



We now turn to our main interest, namely Equation (1.2). We apply the same 



techniques as above in a slightly more complicated setting. After giving all details 
in Smoluchowski's case, we will give a shorter proof and focus on the differences 



with the proof of Theorem 2.2 As above, we can transform the system (1.2) into 



a non-local PDE problem, which we are able to solve, thus obtaining existence and 



uniqueness to (1.2). More precisely, we consider the following system. 



Definition 4.1. Let co(a, m) > 0, a G N, m G N*. We say that a family (cj(a,m)), 



t > 0, a E N, m E N* , is a solution of Smoluchowski's equation (4.1) if 

ft , 



• for every t > 0, (c^, a)^ ds < +oo, 

• for alia eN, m e W and t > 0, 

Ct{a,m) = Co{a,m) + 

' a+l m—1 

+ I ^ a'{a + 2 — a)cs{a, ra)cs{a + 2 — a ,m 

a'=l m'=l 



m 



ds 



(4.1) 



/ aa Cs{a, m)cs{a' , m') ds, 

■^0 a'>lm'>l 



• if (co, a^) < +00, then t h-> (q, a^) is hounded in a right neighborhood of 0. 
Because of the interpretation of variable counting the number of arms a 



particle possesses, it is more natural to state (4.1) in the discrete setting, as in [3]. 
In particular, since at each coagulation two arms are removed from the system, a 
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non-integer initial number of arms would lead to an ill-defined dynamics. One could 
however with no difficulty consider an initial distribution of masses on (0, +00). 
It is easy to see that (q) is a solution to this equation if and only if the function 

+00 +00 

Ux,y) :=5^5^aQ(a,m)x«-ii/"^, (4.2) 

a=l m=l 

defined for t > 0, ?/ G [0, 1] and x G [0, 1), satisfies 



{ks{x, y) - xAs) -Q^i^i y) - Asks{x, y) 



ds, 



(4.3) 



kt{x,y) = ko{x,y) + 
At := kt{l,l) = {ct,a). 

We may solve this PDE with the same techniques as above and obtain the following 
result. 

Theorem 4.2. Consider initial concentrations CQ{a,m) > 0, a G N, m G N* such 
that (co, 1) < +00, Aq := (cq, a) G (0, +00] and with K := (cq, a{a — 1)) G [0, +00] . 
Then K = +00 whenever Aq = +00. Let 

1 



if Ao<K < +00, 
Tgei = s if = +00, 

+00 if K < Aq < +00. 



(4.4) 



When Tgei < +00, this 



Then equation (4.1) has a unique solution defined on K 
solution enjoys the following properties. 

(1) The number of arms At := {ct,a) is continuous, strictly decreasing, and for 
allt>0 



At< 



A, 



l + tAo 



If we set 

then at is given by 



if Aq < +00, 



at = exp 



At < - if Ao = +00. 



(4.5) 



A, ds 



= 1 + Agt for t < Tgel 



and for t >T^ 



gel 



at 



r-i(i + Aorgei + t-Tgei) ifAo< 



-00, 



T-\l + t) 



if Aq = +00, 



(4.6) 



where 

r(x) = 1 + AoTgel + 



dr 



X>1 + AnT 



O-f gel; 



and H : [G'(O), G{1)) 1— ?■ [0, 1) is the right inverse of the increasing function 



G:[0,1)^[G(0),G(1)), Gix):=x--^ 



kn(x, 11 



X G [0, 1) 



(4.7) 
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with G(0) := G(0+) < 0, and 
< G{1) := G(l-) = 



K 



if Aq < +00 
if Aq = +00. 



(2) Let fco be defined as in (4.2), and 



At = {ct,a), 



at = exp 



A, ds 



(3t 



ds. 



a: 



(4. 



Consider 

(ptix, y) := at{x - (3tko{x, y)), 



t>0, x,ye[0,l]. 



Then 



1) attains its maximum at a point it such that 
t < Tge\, d-t = 1; o^nd for t > Tgei, < £j < 1 and 



dx 



(^i,l) = 0. 



In particular, for t > Tgci, d-t is given by 

1 



H 



at 



1. For 



(4.9) 



(4.10) 



where H is the right inverse of the function G defined above. 
• For every y G [0, 1], (l>ti-,y) has a right inverse ht{-,y) : [0, 1] i— )■ [0, 1]. 
(3) The generating function kt defined by (4.2) is given by 



kt{x,y) = —ko{ht{x,y),y) 
at 

for y E [0,1], X E [0, 1] . In particular, for t > 



(4.11) 



atAt = atkt{l, 1) = ko{it, 1), 



At 



koiit,l] 



1 + £ koiis, 1) ds 



(4.12) 



(4) The second moment {ct,a^) is finite on [0,Tgei), infinite on [Tgei,+oo). 



4.1. Proof. The only major difference witli respect to the proof of Theorem 2.2 



IS 



the additional variable y in the generating function kt{x,y). However, the variable 
y plays the role of a parameter in the PDE (4.3), and this allows to adapt all above 
techniques. 

Proof of Theorem \4.S\ The case K < Aq < +oo, for which Tgei = +oo has already 
been treated in [3i Thm. 2], so that we can restrict here to the cases where Tgei < 
+00. When Tgei > 0, Thm. 2 in [3] also shows that at = 1 + A^t on [0,Tgei) (this 
however also requires that {a^,Ct) be bounded in a neighborhood of 0: see point 3 
of the proof of Lemma 2.9). 
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(1) First, by setting ut{x,y) := ath^cptix , y) , y) — ko{x,y), we can see, arguing as 
in points (i)-(ii) of the proof of Lemma 2.9, that for all y G (0, 1] and t > 
there exists it{y) < ^tiy) ^ (0, 1] such that 

atkt{Ux,y),y) = ko{x,y), Vt > 0, y G (0, 1], x e [i'Hy) , lt{y)] (4.13) 

and 4>t{-,y) '■ [^tiy)y^t{y)] ^ [O, l] is a continuous bijection and has a contin- 
uous right inverse ht{-,y) : [0,1] f-^ [i^{y), it{y)]. 




















-L t 1 



Figure 3. 0t(-,l) before and after gelation. The dotted lines repre- 
sent what (pt may look like. The solid one is the actual (pf 

(2) We denote for simplicity 

kt{x) := kt{x,l), := 1), t>0, xG[0, 1]. 

For y = 1, we set £t(l) = £t, i.e. 

l = (j),{et) = at{et-l3tko{it)), t>0. 



Arguing as in points (iv)-(v) of the proof of Lemma 2.9, we can see that 
it = I for all if: < Tgd and £t < 1 for all t > Tg^i- Moreover, t i— )■ is 
continuous and monotone non-increasing. Since (pt is increasing on [0,£t], 
(p'M > 0, i.e. 

1 



so that 



at{it-Ptkom>atG{£t), 



(4.14) 



where we set G{x) := x — x G [0, 1). Notice that 



G'(x) 



1 



{k',{x)f - ko{x)k'^{x) _ ko{x)k'^{x) 



> 0, 



{k',{x)y {k',{x)Y 

since k^ is strictly convex (there is no gelation whenever k'^ = 0). Moreover 
G'(O) < and 

^0 



G(l) 



1 - 



K 



if Aq < +CXO, 



G(l) 



1 if A, 







-l-cxo. 
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Indeed, A;q(1) = K = (cq, a{a — 1)) and, if ko{l) = Aq = +00, then 

ko{x) 



xn k'Jx) 







since, if liminfa;i-i > e > 0, then fco(l) < /i;o(l — 5)e^/^ < +00, for some 
6 > 0, contradicting /i;o(l) = +00. In any case, G has an inverse H, and 
H{l/x) is defined for x G [1 + AoTgd, +00). 
(3) Computing (4.13) at (x, ?/) = 1) we obtain 

d+aj 



ko{(it 

Let us notice that 

<Pt{x) 



atkt{l) = atAt 



X + 



A.. 



dt 



(4.15) 



[X) - 



ds. 



Then by (|4.15|), analogously to (|2.19|) above, 

dt + 0;(£i) d£t = 0;(£i) d£i. 



O = d0i(£t 



koiit 



at 



In particular, for dfj-a.e. t, 0((i't) = 0, i.e. f3t = l/k'Q^it), and therefore 
1 = atiit - (3M^t)) = atG{^t), dit - a.e.t. 



Then, by (4.14), we can write (note that H is well-defined on the considered 
interval) 



< H 



yt > % 



H 



dit — a.e. t. 



Now, by dilsj ), setting A : ]1 + AoTgei, +cx)[^]0, 1[, A{z) := ko [H (i)). 



d+Of 
dt 



< A(aO, Vt > Tgei, 



d+ot 
dt 



A(Q;f), d£j — a.e. t. 



Since at > 1 + AqT^^i for any t > Tgei, we obtain that A;o(^t) < A(aj) < 1 for 
all t > Tgei. In particular, dit is not identically equal to 0. Suppose that for 
some t > Tgei we have (p'tiit) > 0. We set 

s := sup{r < t : 4>'r{ir) = 0} = max{r < t : </)',(£,•) = 0}. 

Then for all r G ]s, t[ we must have (p'riir) > 0. Then for all r G ]s, t[ we have 
ir = is- But, by definition of (3, 



/3r>(3s 



1 



and this is a contradiction. Therefore for all t > Tgei, we have at = A{at) for 
all t > Tgei and the only solution of this equation with a-r ei = 1 + ^o^'gei is 
given by (4.6). 
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(4) In order to prove (4.12), let us note that by the preceding results 

dat 
dt 



atAt = atkt{l, 1) = ko{it, 1), 



A, = ;^loga. = Aiog(l+ / M4,l)d.)=-^f^ 
dt \ Jo J l + £ko{is,l) 



ds 



The rest of the proof follows the same line as that of Theorem |2.2[ □ 

5. The modified version 

Let us finally consider Flory's version of the model with arms. As in the case 
of Flory's equation (3.1), we can consider only initial concentrations cq such that 
^0 = (co,a) < +00. Then, the equation we are interested in is 

^ ^ a+l m—l 

— Ct(a, m) = - a' {a + 2 — a')ct{a' , m!)ct{a + 2 — a' ,m — m!) 

di 2 



a'=l m'=l 



aa'ct{a,m)ct{a' ,m') 



a'>l m'>l 



(5.1) 



^ a\m'>\ / 



act{a, m). 



With the same techniques as above, we can prove the following result. 

Theorem 5.1. Consider initial concentrations Co{a,m) > 0, a G N, m G N* such 
that Aq := (co,a) G (0, +C)o) and with K := {co,a{a — 1)) G [0,+oo]. Let Tgei be 
defined as in (4.4). Then equation (5.1) has a unique solution defined on M^. When 
Tgci < +00, this solution enjoys the following properties. 



(1) We have 



A, 



1 



Mh) 



l + tAo 

where k = 1 for t < Tgei and, for t > Tgei, h is uniquely defined by 
lt = ^^Tirkoik), he [0,1). 



(5.2) 



1+tAo 

Therefore t At is continuous and strictly decreasing on [0, +oo) and ana- 
lytic on ]R"'"\{Tgei}. 

(2) The function (j)t{x,y) = {l + tA())x — tko{x,y) has, for every y G [0,1], aright 
inverse ht{-,y) : [0,1] — )■ [0,lt]- The generating function kt defined in (4.2) 
is given for t > by 

1 



ktix,y) 



-ko{ht{x,y),y). 



1 + tAo 

(3) The second moment {a?,Ct) is finite on ]R+\{Tgei} and infinite at Tgd. 



(5.3) 



Proof. The proof follows the same line of reasoning as the one of Theorem |3.2[ First, 
for every y G [0, 1], 0t(-, y), as defined in the statement, has the following properties: 
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(i) Mo,y)<o, 0t(i,i/)>0f(i,i) = i; 

(ii) For t < Tgei, (f>t{-,y) is increasing , and in particular, there are unique < 

i^y) < h{y) < 1 such that <l)t{^tiy),y) = o and Mh{y),y) = i; 

(iii) For t > Tg^i, 4>t{-,y) is increasing then decreasing for, and in particular, 
there are unique < lf{y) < k{y) < 1 such that = and 

Uh{y),y) = i. 




Figure 4. (pti-, 1) before and after gelation. 
In any case, it is easy to check that for x E [Itiy), h{y)], 

exp ( / Asds) kt{(l)t{x,y),y) = ko{x,y) 







where At is defined by ( 5.2[ ). Then, the properties above show that has a 



right inverse ht defined on [0, 1], and compounding by ht in the previous equation 



shows that (5.3) holds. The other properties then follow easily. □ 

6. Limiting concentrations 

We compute here some explicit formulas for the concentrations and their limit 
for the two models above. In the standard Smoluchowski and Flory cases, particles 
keep coagulating, and they all eventually disappear into the gel: Ct{m) — )■ for every 
m > 1. When the aggregations are limited, there may remain some particles in the 
system, since whenever a particle with no arms is created, it becomes inert, and so 
it will remain in the medium forever. In the following, we consider monodisperse 
initial conditions, i.e. Co(a,m) = fi{a)l^m=i} for a measure fi on N. We also denote 

z/(m) = (m + l)/i(m + 1). 

In [3], it is assumed that u is a probability measure, what we do not require. The 
results of [3] can hence be recovered by taking Ao = 1 below. Now, note the two 
following facts. 

• Equations (4.5) and ( |5.2 ) readily show that 



Coo(.a,m) := hm Ct(a,m) = 0, a > 1, (6.1) 
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that is, only particles with no arms remain in the medium (else, a coagulation 
"should" occur). 

• There is an arbitrary concentration of particles with no arms at time 0, and 
they are the only particles with no arms and mass 1 which will still be in the 
medium in the final state. Hence, the limit concentrations Coo(0, 1) = Co(0, 1) 
have no physical meaning. We will thus only consider Coo(0,m) for m > 2. 

Note now that if at time 0, each particle has zero or more than two arms, then 
obviously, this property still holds for any positive time. Rigorously, this is easy to 



check with the representation formula (4.11) or (5.3). Then, because of (6.1), 

Coc{m) = 

for each m > 2. We thus rule out this trivial case by assuming that 

z/(0) > 0. (6.2) 

This is actually a technical assumption which is needed to apply Lagrange's inversion 
formula in the proof of the following corollaries. We will relate our results to a 
population model known as the Galton- Watson process. For some basics on this 
topic, see e.g. the classic book [2]. The formula providing the total progeny of these 
processes was first obtained by Dwass in f?]. 

6.1. Modified model. 



Corollary 6.1. Let Ct{a,m) be the solution to Flory's equation with arms (5.1) and 
with initial conditions Co(a,m) = fi{a)l^rn=i} with > 0. 

• For allt>0,m>2,a>0, 

, (a + m-2)! t^-^ 

ct(a,m) = — — — -u*"'(a + m- 2). 

^ ' ^ a\m\ (l + tAo)"+™-i ^ ^ 

• In particular, there are limiting concentrations Coo{a,m) = Coo(^)l{a=o} with 

Coo(m) = -^—^u*"'{m - 2). (6.3) 
m[m — 1) 



Proof. With the notation of Theorem |5.1[ we have 



(1 + tAo)ht{x, y) - tyko{ht{x, y)) = x, kt{x, y) = ^ ^^^ yko{ht{x, y)). 
Up to some obvious changes (just replace 1 + t by 1 + tAo), these are precisely 



the equations solved in Section 3.2 of [3] under the assumption (6.2). Theorem 2 
and Corollary 2 therein hence give the desired result (with only 1 + t replaced by 
l + tAo). □ 

If ^0 = 1) which we may always assume up to a time-change, we observe as 
in [3] that 2(m — l)coo(0,m) is the probability for a Galton- Watson process with 
reproduction law started from two ancestors, to have total progeny m. This 
Galton- Watson process is (sub)critical when K := X]a>i '^('^ ~ ^ that is. 



by Theorem 5.1, when there is no gelation, and supercritical when K > 1. Denote 



by py its extinction probability, i.e. the smallest root of k^^x) = x, so p^ = 1 when 
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K < 1 and < 1 when K > 1. Let us compute the mass at infinity, as in [3], by 
writing 

Moo := V mcoo(m) = Coo(l) + V ^-z/*'"(m - 2) 
— ^ — ^ m — 1 

m>l m>2 

Coo(1) + VK«) E ^^z/*'"-'(m-2-a) 
^-^ ^-^ m — 1 

a>0 m>a+2 



n 

a>0 7i>a-|-l 

Now, the Lagrange inversion formula |30] shows that 

n 

is precisely the coefficient of x" in the analytic expansion of around 0, where 
is the unique solution to = Hence 

V - 1 - a) = p^, 

n 

n>a-|-l 

where p,^ is defined above. Note also that Coo(l) = /i(0), so finally 

Moo = Coo(l) + '^(«) ^P^""' = E/^W^*- (6.4) 

a>0 a>0 

The mass at time is Mq = Xl/^l'^); so when there is no gelation, p^ = \ and no 
mass is lost in the gel. When there is gelation, p,y < 1 and the mass Mq — Moo > 
is lost in the gel. By Dwass' formula [7], Moo is also the probability that a Galton- 
Watson process, with reproduction law /i for the ancestor and u for the others, has 
a finite progeny. 

6.2. Non-modified model. 

Corollary 6.2. Let Ct{a,m) be the solution to Smoluchowski's equation with arms 
(4.1) and with initial conditions Co(a, m) = /i(a)l{m=i} with /i(l) > 0. 
• For allt>0,m>2,a>0, 

, (a + m-2)!/3"-i 
alml a? 



where at and (3t are defined in Theorem 4-2 



In particular, there are limiting concentrations Coo(a, "m) = Coo{'m)l[a=o} with 

coo(m) = ^ /3--V-(m - 2) (6.5) 
m[m — 1) 

where /3oo is defined by 
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and c is the unique solution to A;q(c) = ko{c)/c. Moreover, /5oo = 1 when 
there is no gelation, and P^o > 1 otherwise. 



Proof. As for Corollary 6.1, the proof of the formula for Ct{a,m) is the same as in 



[31 Section 3.2], just replacing 1 + tAo by at and t by atPt- So we just have to find 



Now, plj| gives l3t = l//co(^t)) so /St tends to 

1 



the limit of jSt- First (|46]) shows that at +oo, hence, by ( |4.10[ ), = H{0) 



where by definition c := -ff(O) is the unique solution to A;o(c) = kQ{c)/c. Finally, 



when there is gelation, at < I + t after gelation because of (4.6), so by (4.8), 



/3oo > 1. □ 

By a similar computation as above, we may also compute the mass at infinity in 
this case and get 

Moo = $^/u(a) 

a>0 

where c is defined in the corollary. Note that c is the slope of the straight line 
passing by and tangent to the graph of fc, so c > pi,. In particular, less mass is 
lost than in Flory's case. 



A final remark is that despite the striking resemblance between Formulas (6.5) 



and (6.3), the meaning of the factor P^o is unclear. A probabilistic interpretation 



using the configuration model may explain its appearance. 
Acknowledgements We thank Jean Bertoin for useful discussions and advice. 

References 

[1] D. J. Aldous. Deterministic and stochastic models for coalescence (aggregation and coagula- 
tion): a review of the mean-field theory for probabilists. Bernoulli, 5(l):3-48, 1999. 

[2] K. B. Athreya and P. E. Ney. Branching processes. Dover Publications Inc., 2004. 

[3] J. Bertoin. Two solvable systems of coagulation equations with limited aggregations. Ann. 
Inst. H. Poincare Anal. Non Lineaire, 26(6):2073-2089, 2009. 

J. Bertoin and V. Sidoravicius. The structure of typical clusters in large sparse random con- 
figurations. J. Stat. Phys., 135(1):87-105, 2009. 

M. Deaconu and E. Tanre. Smoluchowski's coagulation equation: probabilistic interpretation 
of solutions for constant, additive and multiplicative kernels. Ann. Scuola Norm. Sup. Pisa 
CI. Sci. (4), 29(3):549~579, 2000. 

P. B. Dubovskii and L W. Stewart. Existence, uniqueness and mass conservation for the 
coagulation-fragmentation equation. Math. Methods Appl. Sci., 19(7):571-591, 1996. 
[7] M. Dwass. The total progeny in a branching process and a related random walk. J. Appl. 
Probability, 6:682-686, 1969. 

M. H. Ernst, E. M. Hendriks, and R. M. Ziff. Kinetics of gelation and universality. J. Phys. 
A, 16(10):2293-2320, 1983. 

M. H. Ernst, E. M. Hendriks, and R. M. Ziff. Coagulation processes with a phase transition. 
J. Colloid Interface Sci., 97:266-277, 1984. 
[10] M. Escobedo, P. Laurengot, S. Mischler, and B. Perthame. Gelation and mass conservation in 
coagulation-fragmentation models. J. Differential Equations, 195(1):143-174, 2003. 



POST-GELATION UNIQUENESS FOR COAGULATION EQUATIONS 



31 



M. Escobedo, S. Mischler, and B. Perthame. Gelation in coagulation and fragmentation mod- 
els. Comm. Math. Phys., 231(1):157-188, 2002. 

N. Fournier and P. Laurengot. Well-posedness of Smoluchowski's coagulation equation for a 
class of homogeneous kernels. J. Fund. Anal, 233(2):351-379, 2006. 

N. Fournier and P. Laurengot. Marcus-Lushnikov processes, Smoluchowski's and Flory's mod- 
els. Stochastic Process. AppL, 119(1):167-189, 2009. 

I. Jeon. Existence of gelling solutions for coagulation-fragmentation equations. Comm. Math. 
Phys., 194(3):541-567, 1998. 

N. J. Kokholm. On Smoluchowski's coagulation equation. J. Phys. A, 21(3):839-842, 1988. 
P. Laurengot. Global solutions to the discrete coagulation equations. Mathematika, 46(2) :433- 
442, 1999. 

P. Laurengot. On a class of continuous coagulation-fragmentation equations. J. Differential 
Equations, 167(2):245-274, 2000. 

P. Laurengot and S. Mischler. On coalescence equations and related models. In Degond, P., 
Pareschi, L. and Russo, G. (eds) : Modeling and computational methods for kinetic equations. 
Birkhduser, pages 321-356, 2004. 

F. Leyvraz and H. R. Tschudi. Singularities in the kinetics of coagulation processes. J. Phys. 
A, 14(12):3389-3405, 1981. 

J. B. McLeod. On an infinite set of non-linear differential equations. Quart. J. Math. Oxford 
Ser. (2), 13:119-128, 1962. 

G. Menon and R. L. Pego. Approach to self-similarity in Smoluchowski's coagulation equations. 
Comm. Pure Appl. Math., 57(9):1197-1232, 2004. 

R. Normand. A model for coagulation with mating. Jour. Stat. Phys., 137(2):343~371, 2009. 
J. R. Norris. Smoluchowski's coagulation equation: uniqueness, nonuniqueness and a hydro- 
dynamic limit for the stochastic coalescent. Ann. Appl. Probab., 9(1):78-109, 1999. 
J. R. Norris. Cluster coagulation. Comm. Math. Phys., 209(2) :407-435, 2000. 
G. Stell and R. Ziff. Kinetics of polymer gelation. J. Chem. Phys., 73:3492-3499, 1980. 



R. van der Hofstad. Random graphs and complex networks. Available at http://www.win. 
itue.n l/-rhof stad/NotesRG CN2010.pdf , 

H. J. van Roessel and M. Shirvani. Some results on the coagulation equation. Nonlinear Anal., 
43(5, Ser. A: Theory Methods) :563-573, 2001. 

H. J. van Roessel and M. Shirvani. A formula for the post-gelation mass of a coagulation 
equation with a separable bilinear kernel. Phys. D, 222(l-2):29-36, 2006. 
M. von Smoluchowski. Drei vortrage iiber diffusion, brownsche molckularbewegung und koag- 
ulation von kolloidteilchen. Phys. Z., 17:557-571 and 585-599, 1916. 



H. S. Wilf. Generatingfunctionology. Academic Press, 1994. Also available online at http: 
'/ / www . math . upenn .edu/~wilf/gf ology2 . pdf 

R. M. Ziff. Kinetics of polymerization. J. Statist. Phys., 23(2):241-263, 1980. 

Laboratoire de Probabilites et Modeles Aleatoires, Universite Paris 6 - Pierre 
ET Marie Curie, 175 rue du Chevaleret, 75013 Paris, France 
E-mail address: raoul . normandSupmc . f r 

Laboratoire de Probabilites et Modeles Aleatoires, Universite Paris 6 - Pierre 
ET Marie Curie, 175 rue du Chevaleret, 75013 Paris, France 
E-mail address: lorenzo.zajnbotti@upmc.fr 



